Vibrational properties of the one-component a phase 



S.I.Simdyankin 1 , S.N.Taraskin 2 , M.Dzugutov 1 , and S.R.Elliott 2 

1 Department of Numerical Analysis and Computing Science, 
Royal Institute of Technology, SE-100 44 Stockholm, Sweden. 
2 Department of Chemistry, University of Cambridge, 
Lensfield Road, Cambridge CB2 1EW, UK 
(February 1, 2008) 

A structural model of a one-component er-phase crystal has been constructed by means of molec- 
ular dynamics simulation. The phonon dispersion curves and the vibrational density of states were 
computed for this model. The dependence of the vibrational properties on the thermodynamical 
parameters was investigated. The vibrational density of states of the <r-phase structure is found 
to be similar to that of a one-component glass with icosahedral local order. On the basis of this 
comparison it is concluded that the a phase can be considered to be a good crystalline reference 
structure for this glass. 

(N . 

x> : 

<D . 
P-h; 

\0 ■ 

I. INTRODUCTION 

The local atomic order in disordered condensed materials is well defined and governs many physical propertiesEI. 
Quite often, for a disordered marerial, it is possible to find a corresponding crystal with similar local and even 
^0 , intermediate-range order which give rise to similarities in many structural and dynamical features of these two solids, 
^v^j ' Such a crystal can be regarded as a reference crystalline structure (crystalline counterpart) for the corresponding 
, disordered substance. In some cases, the reference structure can be uniquely defined. The simplest examples are toy 
■ structural models with force-constant and/or mass disorder. In these toy models, the atoms occupy their equilibrium 
positions, at the sites of a crystalline lattice (e.g. simple cubic), which can be considered to be a reference one 
(see e.g.H). Another related example is a binary substitutional alloy, the j-eference system for which is a periodic 
point lattice with one of the two atomic species placed at the lattice sitesB. The disorder in such models does not 
influence the equilibrium positions of the atoms arranged in an ideal crystallineriattice. This makes possible the use 
of approximate analytical approaches (e.g. the coherent potential approximationu) to treat the vibrational properties 
of the models, provided that the vibrational properties of the counterpart crystal are known. 

In amorphous solids, or glasses, the atoms do not occupy the sites of a crystalline lattice, which results in positional 
disorder. For these materials, a choice of a reference structure becomes problematic. Good counterparts can usually 
be found among the crystalline polymorphs having the same (or similar) chemical composition as .the corresponding 
glass. For example, a-cristobalite appears to be a good crystalline counterpart for vitreous silicaQU. 

The, main purpose of this paper is to investigate numerically the vibrational properties of a one-component a- 
phaseoia crystal which-ia-conjectured to be a good crystalline counterpart for a one-component glass with icosahedral 
local order (IC glass)E30. 

The motivation for this choice of a crystalline counterpart of the IC glass is-|the following. The computational 
model of the IC glass is based on a simple empirical pair interatomic potential!!] resembling the effective interionic 
potentials conjectured for simple metallic glass-forming alloysO. The use of the same potential allows_us to construct 
models of bec and tx-phase crystals that are stable for a wide range of thermodynamical parametersEj. Of these two 
crystalline structures, the a phase is expected to be a good reference structure for the IC glass because of the following 
reasons: The supercooled IC liquid (where the interactions between atoms are,-described by the same potentiaO) 
undergoes a transition either to the IC glass or to a dodecagonal quasicrystalli-3 depending on the quench ratetJ. 
This quasicrystal has similar local structural properties with the IC glassliil. However, the absence of global periodic 
order in the quasicrystalline phase makes the analysis of its vibrational properties a task of comparable complexity to 
that for thepglass itself. The a phase is one of the closest low-order crystalline approximantsE3 for this dodecagonal 
quasicrystalclj, which means that these two (crystalline and quasicrystalline) structures are built up from the sa: 
structural units. This implies that the IC glass and the a phase, being both tetrahedrally closed-packed structure: 
are nearly isomorphous in terms of local order. 

Knowledge of the vibrational properties of the er-phase crystal allows for a direct comparison with those of the IC 
glass. The apparent similarity in the vibrational densities of states of these two structures gives stronger support for 
the choice of this crystalline counterpart for the IC glass. 
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The cr-phase structure, used in our computations, has been obtained by means of molecular dynamics simulation 
with the use of an interatomic pair potentiate. The vibrational properties have been investigated by using both a 
the normal-mode analysis and by computing the spectra of appropriate time-correlation functions. 

The paper is arranged as follows: The cr-phase structure is described in Sec. O. The model and technical details 
of the calculations are presented in Sec. III. In Sec. |y| we present the results of the simulations. Some concluding 
remarks are contained in Sec. 0. 



II. THE a PHASE 



In this Section, we review the known structural and dynamical properties of the a phase. 



A. Structure 

The a phase beloii&s to an important class of tetrahedrally close-packed crystallographic structuresi, viz. the 
Frank-Kasper phasea3 : El. The first coordination shells of the constituent atoms in these structures form triangulated 
(Frank-Kasper) polyhedra composed entirely of slightly distorted tetrahedra. The four possible coordination numbers 
(Z) in these structures are Z = 12, 14, 15 and 16. The least distorted tetrahedra are found in icosahedra (Z12 
polyhedra). Structures of small clusters of atoms interacting via pairwise-.central potentials favor icosahedral ordercl 
as having the lowest energy. The prototype a phase structures are /3-tO and Cr4sFe52t3. There are 30 atoms per 
tetragonal unit cell (tP30) with c/a w 0.52, where c and a are the dimensions of the cell (lattice parameters). The 
space group of this phase is P^/mnm. There are 10 atomSj-with the coordination number 12, Z12 or icosahedra, 
16 Z14 atoms and-,4 Z15 atoms. The —72° disclination linesEl form a network (a major skeleton, in the parlance of 
Frank and Kasperuu), where rows of 214 atoms parallel to the tetragonal c-axis thread planar networks of Z14 and 
Z15 atoms. A projection of the cr-phase structure down the c-axis is shown in Fig. [|. 

The-Fxank-Kasper phases share their significant geometrical property of icosahedral local order with simple metallic 
glassesBtj. Some liquid alloys which form Frank-Kasper phaseSrhave a tendency to freeze into metastable amorphous 
structures (metallic glasses) when quenched sufficiently rapidlylia. It is now generally well accepted that, at least in 
the case of metallic alloys of simple constitution, glass formation is caused by the incompatibility of local icosahedral 
coordination with the translational symmetry in Euclidean space feeometrical frustration)!! There exist statistical 
mechanical arguments in favor of this scenario of glass formationO based on a Landau free-energy analysis. The 
average coordination numher, Z, in the Frank-Kasper phases (13.333 < Z < 13.5) is very close to that of a sphere- 
packed "ideal glass" modelEll = 13.4). In a sense, this "ideal glass" could be regarded as a Frank-Kasper phase 
with an infinitely large unit celloEJ. Thus the class of Frank-Kasper phases is a natural choice for reference crystalline 
structures for metallic glasses of simple constitution. 

From a structural point-, of view, the a phase, can be also regarded as a crystalline low-order approximant for 
dodecagonal quasicrystalscj. Such quasicrystalsli3, morphologically close to Frank-Kasper phases, represent an alter- 
native class of noncrystallographic structures which combine icosahedral local order with non-translational long-range 
order manifested by infinitely sharp diffraction peaks. 



B. Dynamics 

Similarities in the local structure of metallic glasses and Frank-Kasper phases are reflected in the dynamical prop- 
erties of these materials. 

The available data about the vibratinnal, dynamics of Frank-Kasper phases is limited to some of the Laves phases, 
a subclass of the Fp,»k-Kasper phasestaO. For instance, the, similarity between the phoitopidispersion relations of 
the Mg7oZn3o glasse3£3 and those of the Laves phase MgZn^ was emphasized by HafnerLxEa. 

Another interesting aspect of the dynamics of Frank-Kasper phases is related to the appearance of soft vibrational 
modes in these materials. Such a soft low-frequency optic mode at the T-point (the origin of the reciprocal lattice) 
has been found numerically in the same Laves phase MgZn2E.ll. The frequency of this mode decreases with increas- 
ing pressuxe Accompanied by volume. compression) and eventually becomes negative, indicating a structural phase 
transitionE3~EJ. The authors of Ref.Efl applied a group-theoretical analysis and demonstrated that the polarization 
vector of the soft optic mode in MgZn 2 is determined by the structure symmetry and is independent of interatomic 
interactions. This suggests that the soft-mode character of some vibrations is a generic property of Frank-Kasper 
phases. However, no soft-mode behavior was observed in the isostructural CaMg2u where the mass ratio of the 
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constituent elements is reversed with respect to MgZn2. Hafnerlll suggested that the soft modes in MgZn2 should 
be attributed to the relatively large mass of the Zn atoms. This is an example of the chemical composition of the 
materials introducing considerable difficulties into the analysis of the interplay of the structure and the dynamics. 

A numerical simulation of a one-component Frank-Kasper phase allows us to eliminate this uncertainty. We in- 
vesti gate th e behavior of the lowest-frequency optic modes in the one-component a phase with variable pressure in 
Sec. 
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III. METHODS OF COMPUTATION 

We have constructed a thermodynamically stable structural model of the one-component a phase by means of 
classical molecular dynamics simulation. In this case, success strongly depends on the choice of the interato: 
potential for the modeL-|For example, the Lennard-Jones potential, widely used in creating simple models of liquid: 
glassesEH, and crystals^, is not suitable for this purpose, because the a phase is not stable for it in the range of 
thermodynamical parameters investigated below (the stable phase is the fee lattice). Instead, we use a pair interatomic 
potential suggested in Ref.E3 and show that-it is possible to construct a one-component a phase which is stable for a 
wide range of thermodynamical parameters!!! 

A. Model 

As a mathematical model for the study of atomic dynamics in the crystalline cr-phase structure, we consider a 
classical system of N identical particles-interacting via a spherically symmetric pair potential. 

The pair potential used in this studyEUl is designed to favor icosahedral local order. The main repulsive part of this 
potential is identical to that of the Lennard-Jones potential ULj(r) fn, 4e[(er/r) 12 — (er/r) 6 ]; therefore all the quantities 
in these simulations are expressed in reduced Lennard-Jones unitso, i.e. with tr, e, and tq = (ma 2 /e) 1 / 2 chosen as 
length, energy, and time units. To convert the reduced units to physical units one can refer to argon (m = 39.948 
a.m.u.) by choosing the Lennard-Jones parameters a = 3.4 A and e/fcs = 120 K. In this case, oik frequency unit 
i/q = Tq 1 corresponds to 0.4648 THz. The analytical expression defining the potential is given in Ref.EJl This potential 
resembles those for simple glass-forming metallic alloyst3 with only the first of the Friedel oscillations being retained 
(seeFig.|). 

We use conventional molecular dynamics simulations^! in which the Newtonian equations of motion are solved using 
a finite-difference algorithm with time step equal to 0.01 while the particles are enclosed in a simulation box of volume 
V with periodic boundary conditions. In this case, the total energy E is a constant of motion and time averages 
obtained in the course of simulations approximate the ensemble averages in the microcanonical (constant-AX^i?) 
statistical ensemble. Wavevectors, 

Q = n x Qx,o + n y Q y fi + n z Q z ,o (3-1) 

with n x , n y and n z being integers, compatible with periodic boundary conditions are multiples of the three fundamental 
wavevectors Q x ,o = ^(1,0,0), Q y! o = ^"(0,1,0) and Q 2i0 = 1^(0,0,1), where L x , L y , and L z are the dimensions 
of the (tetragonal) simulation box. In order to have a sufficient number of allowed wavevectors within the Brillouin 
zone, the sample dimensions must be sufficiently large. The time-correlation functions resulted from the molecular 
dynamics simulation reported below were obtained for a system of 20580 particles (7 x 7 rX.14 unit cells with 30 
atoms per unit cell). This relatively large system size also gives sufficient statisti£aLaccuracyE3. Where necessary, we 
performed simulations in other ensembles by modifying the equations of motiontiZrLj. 

The arrangement of atoms in a unit cell of the model cr-phase structure used in our computer simulations is shown 
in fig. 0(b). The optimal (with respect to an energy minimization) c/a ratio was taken to be equal to 0.5273. To 
determine the optimal model structure, it was sufficient to perform molecular dynamics simulations with the number 
of particles in the sy stem equal to N = 1620 (3x3x6 unit cells) . Details of the preparation of this atomic configuration 



are given in section IV A 



B. Time-correlation functions 



A straightforward method to analyze the vibrational dynamics in a molecular dynamics model is to imitate inelastic 
neutron scattering experiments by calculating the dynamical structure factor S^Q, w)Ej, proportional to the neutron 
scattering cross-sectionEZtES, which is the spectrum of the density autocorrelation function: 
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F(Q,t)=<p(Q,t)p(-Q,0)> (3.2) 



where 

N 



p(Q,t)=^exp(-iQ-r fe (0) (3.3) 



k=l 



is the Fourier transform of the local particle densityil, N is the number of particles in the system, Tk(t) is the position 
vector of particle k, and the wavevector Q takes on the values according to Eq. (3.1). A longitudinal phonon is 
associated with a maximum in S(Q,u>) at a fixed Q. In order to get-uiformation about the transverse modes from 
5(Q,w), one has to select wavevectors outside the first Brillouin zoneSa. .— . 



In a more convenient way, the vibrational modes can be studied using the current autocorrelation functio: 



Ce(Q,*) = <Je(Q,*)je(-Q,0) > (3.4) 



where 

N 

je(Q J t)=E(e-Vfc(t))exp[-*Q.r fc (i)] (3.5) 
fe=i 

is the Fourier transform of the local current, Vfc(t) is the vel ocity of particle k, and e i s the unit polarization vector. 

Note that for the longitudinal polarization, e || Q, Eq. ( |3.4[ ) can be obtained from d3.2| ) by double differentiation with 
respect to time. For the transverse-current correlation functions, the polarization vector must be chosen consistently 
with the lattice symmetry. 

At a temperature T, the vibrational density, of states g(ui) can be calculated as the Fourier transform of the 
normalized velocity autocorrelation functionc3H: 

We computed the time-correlation functions using the overlapped data collection technique!!. The number of 
overlapped measurements used for statistical averaging was about 10000. The time origins of the measurements were 
separated by 0.2 r.u. (20 time steps). 

In order to reduce the finite-time truncation effects in the spectra of the time-correlation functions, we used a 
Gaussian window function with the half-width equal to 3 r.u. 



C. Normal- mode analysis 

To calculate the dispersion relations in the harmonic approximation, we used the standard method based on 
diagonalization of the Fourier transformed dynamical matrixcil. From the known dispersion relations, w s (Q), s — 
1, 2, . . . , 90, the vibrational density of states can be computed by integration over the first Brillouin zone according to 

g( LJ ) = ——^J2 S(u-cj a {Q))dQ (3.7) 
r \ Z7T ) s= iJbz 

where the sum is over all r dispersion branches and v stands for the volume of the unit cell. In the computations, 
a Gaussian function with small but finite width is substituted for the ^-function. The half-width of the Gaussian 



function in our computations was equal to about 0.05. Formally, the wavevector Q in Eq. (3.7) is a continuous 
variable, but in the simulations the integral was estimated by a sum over a uniform rectangular grid of 100 x 100 x 100 
points in the first Brillouin zone. 
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IV. RESULTS 



A. Optimization of the structure in the a phase 

In this subsection, we describe the method of construction of the thermodynamically stable model of the a phase 
and analyse the range of its stability. 

There are several ways to obtain numerical values of the atomic coordinates in the a phase. One way is to 
use those available for Cr4sFe52c3. Alternatively, a unit cell of the cr-phasc structure can be constructed either by 
manipulating the kagome tiling according to the algorithm given by Frank stH|d-iKasper in refE or by stacking the 
square and triangular basic elements of the dodecagonal quasicrystal modelE3cj into the 3 2 ,4,3,4 square-triangle 
netB (see fig. 0(a)). The arrangements of the atoms resulting from these constructions do not correspond exactly, 
although the difference is rather small - the root-mean-square distance between the corresponding atoms in different 
configurations is of the order of a few percent of the root-mean-square distance between different atoms in the same 
configuration. In either case, the resulting structure is an approximate one in the sense that the atomic positions 
do not correspond to a minimal potential-energy configuration for a given interaction potential. To obtain the true 
structure corresponding to the potential, the approximate configuration must be relaxed by a molecular dynamics 
program. Moreover the c/a ratio is slightly different for different natural tr-phase crystals, which means that this ratio 
is not uniquely defined. 

The atomic configuration of the cr phase used in this study was prepared as follows. A sample of N = 1620 
particles (3x3x6 unit cells) was constructed by filling a tetragonal bcp^pf, appropriate dimensions with cr-phase unit 
cells. We used the unit-cell atomic configuration suggested by Gahlerli3e3 with c/a — 0.5176. The number density, 
p = N/V, V — L x L y L z , of this atomic configuration was p = 1.0048. This configuration was then used to provide 
the initial atomic coordinates for variable-shape NST (constant number of atoms, pressure tensor, and temperature) 
molecular dynamicsES, the run performed at S = and T = 0. This procedure is equivalent to a potential-energy 
minimization by the steepest descent method under the condition of independent pressure balance in each spatial 
dimension. The variable-shape NST run resulted in an ideal crystalline structure for which the fractional coordinates 
of the atoms in all unit cells were identical within the precision of the calculations. The structure of the a phase thus 
obtained is characterized by the minimum potential energy per atom, J7 m ; n = —2.5899, with respect to variations of 
thermodynamical parameters. In order to check this, we have performed similar NST runs at different pressures and 
indeed found that the energy is minimum at zero pressure (see Fig. |3^,,b) . The density for the optimal structure has 
been found to be.p = 0.8771 and the ratio of the lattice parameters is c/a = 0.5273, close to that of /3-U (c/a = 0.5257, 
at T — 720°C)c3. The potential-energy minimum for the bec structure at the same density with the same potential 
was U m [ n = —2.6148. At zero pressure, the density of the bec structure is p = 0.8604 and the minimuip-|potential 
energy per atom is U m i n = —2.6357, i.e. in both cases lower than for the a phase. It was shown earlierll3 that the 
potential energy per atom for the cr phase becomes lower than that for the bec structure at the same density as 
the temperature increases. This is consistent with the fact that natural crystalline cr phases are stable only at high 
temperatures. They undergo a solid-solid phase transition to a simpler crystalline phase as the temperature decreases. 
In our simulations, however, the system was stable in the range of temperatures T < 0.9 for as long as t lun = 5000. 

We have also investigated the thermodynamical stability of the cr phase under variable pressure. We have found that 
the cr phase is stable for pressures in the range —5 < P < 12. At high pressures P > 12, a structural transformation 
occurs, resulting in the fee structure. The phase diagram of the IC potential is not known at preseat. We can 
expect that at densities greater than the triple-point density for the Lennard-Jones system, p « 0.85E3, the solid- 
fluid coexistence curve for the IC system is close to that for the LJ system. We can only estimate that the melting 
temperature at p = 0.8771 is about 0.8 < T < 0.9 from the fact that the cr-phase crystalline structure is stable at 
T = 0.8 on the time scale of our computations. No diffusion was observed at temperatures up to T = 0.8. At T = 0.9 
the system stayed in a metastable superheated state for about t TUa = 5000, after which it melted. 

We used the density and the c/a ratio obtained from the NST run to perform NVE (constant number of atoms, 
volume, and total energy) molecular dynamics runs starting from the three configurations mentioned above and scaling 
the velocities to zero at each time step, which is also equivalent to a potential-energy minimization by the steepest- 
descent method. The same was done for one instantaneous configuration corresponding to the temperature T = 0.8. 
The configurations resulting from this procedure were identical, which is an indication that there is a well-defined 
potential-energy minimum corresponding to a unique crystallographic arrangement of atoms within the cr-phase unit 
cell. This structure, scaled so that a = c = 1, is shown in fig. |l](b). The atomic layers with z close_to 0.25 and 0.7IL 
which are not closely packed, show a small but significant puckering - an effect present in the /3-UL3 and Cr4gFe520 
structures. 
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B. Vibrational dynamics 



Above, we have discussed the similarities in the local structure of the a phase and the IC glass. These similarities are 
expected to cause the vibrational spectra in these two materials also to be similar. In order to check this assumption, 
in this subsection we investigate the vibrational properties of the a phase and compare the vibrational spectra of this 
crystal and the IC glass. 

If the vibrational spectra are similar, the a phase can be considered to be a good crystalline reference structure for 
the IC glass. One consequence of these similarities is that we can use the data about the vibrational properticies of 
the a phase crystal to explain the nature of the vibrational excitations in the corresponding amourphous structure. 



1. Phonon dispersion in the a phase 



The a phase has 30 atoms per unit cell which result in 3 acoustic and 87 optic branches, as shown in Fig. The 
vibrational density of states (VDOS) obtained by integration over the first Brillouin zone (see Eq. (j|)) is shown on the 
right-hand side of Fig. ||. The linear dispersion of the acoustic branches in the low-frequency range (u> < 4) results in 
the Debye law for the VDOS, g(u>) = 3w 2 /wj3, with the Debye frequency equal to lo-q ~ 23.89. The Debye frequency 
has been estimated from a fit of the initial part of the VDOS by a parabolic function. The optic branches are densely 
distributed above the acoustic part. There are no large gaps in the spectrum, which is a consequence of tight binding 
and the mutual penetration of the basic structural units (Frank-Kasper pelyhedra) in the a phase. |-J.n other words, 
there are no isolated structural units, like molecules in molecular crystals^ and crystalline fullerensc3, or tetrahedra 
in silicafl, the vibrations of which form separate optic bands. At some of the zone boundaries (e.g. the X point; see 
Fig. ^) the dispersion curves do not show zero derivatives. This is because the space group PA2 /.amm of the a phase 
is nonsymmorphic, i.e. it contains nonpoint symmetry elements involving fractional translationsoo. 



2. Comparison with the IC glass 



An informative characteristic of the vibrational dynamics in the IC glass which can be compared with t he a phase 
is its VDOScZl (see Fig. ||) which can be easily obtained from the velocity autocorrelation function Eq. (3.6). The 
VDOS for the bec lattice is also presented for comparison in the figure (the dashed line). We can clearly see that the 
frequency range of the whole spectrum is the same for the a phase and the IC glass but differs for the bec lattice. The 
shape of the IC-glass spectrum mainly reproduces the basic features of the a phase spectrum and can be imagined as 
a superposition of broadened (by disorder) crystalline peaks. This is a consequence of the presence of a large number 
of optic modes in the vibrational spectrum of the sigma-phase structure located in the same frequency region as the 
whole spectrum of the IC glass. Therefore, the similarities in the VDOS of the a phase and the IC glass strongly 
support the choice of the a phase as a crystalline counterpart. 



3. Soft modes in the a phase 



As was mentioned in Sec. II B , an interesting feature of atomic dynamics in the Frank-Kasper phases is related 
to the appearance of low- frequency soft modes. We have investigated whether a soft mode appears in our model of 
the cT-phasc structure. For this purpose, we followed the evolution of the vibrational spectrum with variable pressure 
(see Fig. ^ and, indeed, found that one of the lowest-frequency optic modes (doubly degenerate) in the T-point 
shows soft-mode behavior. The frequency of this mode decreases both with decreasing and increasing pressure (see 
Fig. ^) . The decrease of the mode frequency at negative pressures is not surprising and reflects the softening of the 
whole vibrational spectrum (see Fig. Ha). However, with increasing pressure, the whole spectrum is shifted to higher 
frequencies (see Fig. ^c), while the frequency of the soft mode (a small peak around u> = 3.5 in Fig. |^c) moves in the 
opposite direction, approaching zero and thus indicating a structural instability (structural phase transition to the 
fee lattice) at a critical pressure P* ~ 12.5 (see Fig. ^). Around this value of the pressure the structure of the a phase 
becomes extremely unstable and an investigation of the details of atomic motion requires a thorough analysis. We 
hope to address this point in another study. 
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4- Anharmonicity in the a phase 



One of the interesting questions concerning the vibrational dynamics of the a phase is related to the range of appli- 
cability of the harmonic approximation for the lattice vibrations. We are able to anwer this question by investigating 
the vibrational spectrum using the velocity autocorrelation function with increasing temperature and comparing it 
with the results of the normal- mode analysis (harmonic approximation) . 

To assess the degree of temperature-induced anharmonicity, we computed the dispersion relations for the symmetry 
direction [001] (Q || c, TZ in Fig. |J) at different temperatures by using both these techniques (molecular dynamics 
and normal-mode analysis) and compared the results. These are shown in Fig. || for several low- and high-frequency 
dispersion branches for two temperatures T = 0.01 and T = 0.8. At intermediate frequencies, the density of dispersion 
branches is so high that a comparison between the results of the two methods of calculation of dispersion relations is 
hardly possible, m ainly because of the finite width of the respective peaks in the spectra of the current autocorrelation 
functions (see sec. [ill B| ) . Due to the fact that the cr-phase space group P^/mnm is nonsymmorphic, i.e. it contains 
nonpoint symmetry elements involving fractional translationsCJ, the phonoxtdispersion relations, derived from the 
peak positions in C;(Q, to) and Ct(Q, u>), appear in the extended zone schemeBj. The optic modes cannot be measured 
in the vicinity of the origin of the first Brillouin zone Q — (Q — |Q|), because this long- wavelength limit corresponds 
to motion of the system as a whole, forbidden by the periodic boundary conditions. Information about these modes is 
available at the boundaries (Q = 2ir/c, 6tt/c) and at the origin (Q = 4ir/c) of the second extended Brillouin zone. The 
molecular dynamics results for the dispersion relations at T = 0.01 are adapted from the second extended Brillouin 
zone. To make possible the comparison with the results obtained in the harmonic approximation, the data in the 
region Q s [5tt/c, 6tt/c] were folded with respect to Q — 5ir/c into the region Q £ [4tt/c, Stt/c], which corresponds to 
half of an irreducible zone. From these results we can see that the harmonic approximation works quite well at the 
low temperature of T = 0.01. At T = 0.8, only the acoustic branches could be resolved without ambiguity. Therefore, 
for this temperature, we used the data available from the first Brillouin zone. 

One important signature of temperature-induced anharmonicity is a softening of the acoustic modes, i.e. a lowering 
of the acoustic branches with respect to those calculated in the harmonic approximation which occurs as the tem- 
perature increasesEj. This effect can be clearly seen for T = 0.8 in Fig. |8|. In accordance with this observation, the 
vibrational density of states for this temperature, shown in Fig. 0, exhibits the presence of excess states with respect 
to the harmonic approximation. A deviation from the harmonic approximation in g(u>) at low frequencies (see the 
inset in Fig. [)]) starts to be noticeable at a temperature of about T = 0.2. Therefore, we can conclude the the lattice 
dynamics in the a phase is harmonic in a wide range of temperatures T < 0.2. 

Finally, we would like to note the similarity between the high-temperature VDOS for the a phase and the low- 
temperature VDOS for the IC glass (see Fig. |Io|). Since g(uS) for the glass is only slightly temperature dependent, 
we show it only for T — 0.01. The fact that the densities of states for the glass and the high-temperature a phase 
in fig. [l0| are remarkably similar clearly indicates that the effect of the thermally-induced dynamical disorder in the 
crystalline structure on the vibrational spectrum is similar to that of the configurational disorder characteristic of the 
amorphous structure. 

V. CONCLUSIONS 

In this paper we have studied the structural and vibrational properties of a a- phase crystal. First, we have shown 
that it is possible to construct a structural model of a one-component a phase by means of molecular dynamics 
simulations using an appropriate pair potential. This cr-phase structure is stable in a wide range of thermodynamical 
parameters. Our model of the a phase contains only one atomic component. This is important in understanding the 
role of topological icosahedral order alone on the structural and dynamical properties and avoids the effects arising 
from the presence of different atomic species. 

Second, we have investigated atomic vibrational dynamics of the a phase. In particular, we have found the range 
of applicability of the harmonic approximation in a description of atomic dynamics. We have also demonstrated the 
existence of soft modes in the a phase which leads to a structural phase transformation with increasing pressure. 

Third, we have demonstrated that the a phase is a good crystalline counterpart of the IC glass. This has been done 
on the basis of a comparative analysis of the vibrational dynamics (vibrational density of states) . 

We think that the results on the vibrational properties of the er phase (discussed above can be used in an analysis 
of the peculiar vibrational properties of the IC glass (e.g. the Boson peaku). We also believe that the computational 
data of the vibrational properties of the a phase could be of value for metallurgy where this phase has received much 
detailed attention, chiefly ,-hecause of the detrimental effect which the formation of this phase has on mechanical 
properties of certain steelsEl 
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FIG. 1. Projection down the c-axis of the cr-phase structure: (a) 3 2 , 4, 3, 4 net (the numerical symbols are Schlafli symbolffl, 
specifying the number and sequence of various polygons around each vertex) . The dashed square outlines a unit cell; (b) Atomic 
arrangement in one unit (cubic) cell. 0: z = 0, ♦: z = 0.5, •: z = 0.2499. x: z = 0.2501, O: z = 0.7499, □: z = 0.7501, O 
Z12 atoms, Q: Z15 atoms. The rest of atoms are Z14. Multiplying a by the proper c/a ratio gives a tetragonal unit cell with 
a = l. 
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FIG. 3. (a) Minimal potential energy of the a-phase structure as a function of density, (b) Pressure as a function of density. 
The dots show the data points and the solid curves are obtained by a cubic interpolation. 
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FIG. 4. Phonon-dispersion relations along the symmetry directions (a) and vibrational density of states (b) evaluated by 
a normal-mode analysis in the harmonic approximation. The symmetry points on the surface of the first Brillouin zone are 

r = (0,0,0), m = (j, i,o), r = (i 5, i), z = (o,o, f ), a = (j,o, i), x = (5,0,0). 
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FIG. 8. Some phonon-dispersion curves calculated by a normal-mode analysis in the harmonic approximation (solid lines) 
and derived from the spectra of longitudinal and transverse particle-current autocorrelation functions (symbols) for the cr-phase 
structure. O: longitudinal phonons; □: transverse phonons at T = 0.01. •: longitudinal phonons; ■: transverse phonons 
at T = 0.8. The linear size of the symbols is approximately equal to the width of the spectral peaks. The direction of the 
wavevector Q = [0, 0, £], FZ, is parallel to the c-axis. 
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FIG. 9. Vibrational densities of states calculated by normal-mode analysis (solid line) and from the velocity autocorrelation 
function at different temperatures (broken lines) for the cr-phase structure. Dotted line: T = 0.01, dashed line: T = 0.2, 
dashed-dotted line: T = 0.8. The inset shows the low-frequency part of the spectrum 
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FIG. 10. Vibrational densities of states calculated from the velocity autocorrelation function for the <r-phase structure and 
the corresponding glass at the reduced temperatures shown. Both results are obtained at the same density p — 0.8771. 
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